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abstract  -  ECG  lead  selection  and  reconstruc¬ 
tion  were  investigated  in  the  present  study  us¬ 
ing  ECG  source-to-measurement  transfer  matrices 
computed  in  inhomogeneous  and  homogeneous  con¬ 
ductor  thorax-heart  models,  which  represent  end- 
systolic  and  end-diastolic  cardiac  phases.  ECG 
leads  were  selected  from  a  set  of  120  leads  for  body 
surface  potential  mapping.  The  transfer  matrices 
of  the  120-lead  system  were  reconstructed  from  a 
transfer  matrix  of  selected  leads.  Effects  of  chang¬ 
ing  the  cardiac  phase  and  the  conductivities  on  re¬ 
construction  performance  are  reported  here.  The 
number  of  selected  leads  well  matched  to  the  results 
of  earlier  works  conducted  by  other  researchers. 
Keywords  -  Body  surface  potential  mapping,  ECG 
lead,  cardiac  phase,  volume  conductor  model 

I.  Introduction 

Body  surface  potential  mapping  (BSPM  or  BSM)  is 
a  method  to  display  electric  field  potential  distribu¬ 
tion  on  the  torso  surface  generated  by  the  cardiac  ac¬ 
tivity.  From  generated  maps,  features  related  to  di¬ 
agnostic  categories  are  observed  [1].  There  are  BSM 
measurement  configurations  mostly  varied  by  the  num¬ 
ber  of  ECG  leads  and  electrode  positions,  for  example 
[2,  3,  4]  and  more  summarized  by  [5].  While  there  is 
not  a  standard  configuration  for  BSM,  several  compu¬ 
tational  studies  have  been  reported  aiming  at  speci¬ 
fying  necessary  electrode  positions  or  at  reducing  the 
number  of  positions.  Some  of  them  are,  say,  dataset 
approach,  e.g.,  Barr  et.  al.  [6],  Lux  et.  al.  [7],  and  Ko- 
rnreich  et.  al.  [8],  which  use  real  ECG  data,  and  the 
others  are  model  approach,  e.g.,  Dossel  et.  al.  [9],  which 
use  ECG  source-to-measurement  transfer  system  com¬ 
puted  in  the  thorax-heart  models.  The  present  study 
is  a  model  approach  and  applies  a  lead  selection  al¬ 
gorithm,  which  is  mathematically  the  same  as  that  in 
Lux  et.  al.’s  dataset  approach,  to  a  matrix  of  source- 
to-measurement  transfer  coefficients.  The  ECG  source 
in  Dossel  et.  al.’s  model  was  the  epicardial  potentials, 
while  in  the  present  study’s  are  dipole  current  sources 
distributed  in  the  entire  heart. 

Since  the  potentials  on  BSM  are  detected  by  unipo¬ 


lar  leads,  selection  of  electrode  positions  and  that  of 
ECG  leads  can  be  considered  as  questions  equivalent 
to  each  other  if  we  are  allowed  not  to  be  concerned 
with  the  reference  electrode  configuration. 

A  uniqueness  of  the  present  study  may  be  that  lead 
selections  and  reconstruction  coefficients  are  computed 
in  models  of  different  cardiac  phases  and  with  differ¬ 
ent  conductivity  modelings,  and  those  obtained  in  each 
model  are  also  applied  to  the  other  models. 

II.  Study  Material 

Four  thorax  models  with  abbreviations  used  in  the 
present  study  are  listed  below. 

•  i-sys,  Inhomogeneous  conductor  model,  end- 
systolic  cardiac  phase, 

•  i-dia,  Inhomogeneous,  end-diastolic, 

•  h-sys,  Homogeneous,  end-systolic, 

•  h-dia,  Homogeneous,  end-diastolic. 

They  were  inhomogeneous  models  and  infinite  homo¬ 
geneous  models.  With  each  conductivity  modeling, 
two  thorax-heart  models  were  constructed;  one  had  the 
heart  of  the  end-systolic  cardiac  phase,  and  the  other 
had  the  heart  replaced  with  that  of  the  end-diastolic 
phase.  In  each  of  the  inhomogeneous  models,  ECG 
source-to-measurement  transfer  matrix,  C,  was  com¬ 
puted  using  the  finite  difference  method  (FDM).  C 
transforms  dipole  current  source  vectors  at  all  FDM 
node  points  forming  the  heart  model.  The  numbers  of 
such  nodes  were  12032  and  13120  in  the  end-systolic 
heart  model  and  in  the  end-diastolic  model,  respec¬ 
tively.  The  full-BSM  lead  system  had  120  unipolar 
leads  with  the  Wilson  central  terminal  as  the  reference 
level  including  117  body  surface  leads  and  three  limb 
leads  [2].  Therefore,  the  matrix  size  of  C  was  120-by- 
M  where  M  was  12032  x  3  in  the  systolic  phase  model 
and  13120  x  3  in  the  diastolic  phase  model.  A  vector  u 
consisting  of  source  coefficients  of  all  dipole  sources  is 
transformed  to  a  120-length  vector  v  consisting  of  120 
lead  voltages: 

v  =  Cu  (1) 
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Each  row  in  C  transforms  u  into  an  ECG  voltage  of 
one  of  the  BSM  leads. 

The  transfer  matrices  in  the  homogeneous  models 
were  computed  using  the  dipole  field  equation  [10],  but 
the  source  points  and  the  electrode  locations  were  the 
same  as  in  the  inhomogeneous  models. 

III.  Method 

The  lead  selection  and  reconstruction  method  here 
is  based  on  the  Gram-Schmidt  Orthogonalization 
(GSO)  [11],  and  it  is  an  iterative  algorithm  ap¬ 
plied  to  /  vector  entries  a,— i in  a  matrix  A  = 
[ai  &2  ■  ■  ■  a/]  which  is  the  transpose  of  an  ECG 

source-to-measurement  transfer  matrix,  i.e. ,  A  <—  C*. 
/  is  the  maximum  number  of  ECG  leads;  namely, 
/  =  120.  A  can  be  given  as  a  matrix  product 

A  =  QR  (2) 

where  the  columns,  q,— i of  matrix  Q  are  mutually 
orthogonal  basis  vectors  of  A,  and  R  is  an  I-by-I  ma¬ 
trix.  Let  O  denote  the  inverse  matrix  of  R  so  that  Q 
is  given  as 

AO  =  Q  (3) 


where  t'N.i  is  the  element  on  IVth  row  in  ith  column  of 
R,  and  q^e-^  is  the  inner  product  of  qjv  and  e^. 

As  seen  in  (7),  q„=i:jv  are  copies  of  ei"-11  which  are 
transformations  of  N  columns  in  A  whose  indices  are 
first  N  entries  in  the  selection  sequence  S. 

0(N)  consisting  of  N  columns  left-most  in  O  gives 
<ln=l:N  (see  (3)): 

AO(JV)  =  Q(JV)  (9) 

Using  the  corresponding  N  rows  in  R  as  the  entries  in 
a  matrix  RW,  the  reconstruction  matrix,  denoted  by 
HW,  is  defined  as  =  O^R^.  trans¬ 

forms  a.;,  whose  indices  i  are  first  N  elements  in  S,  to 
A(A  ):  AH,Ar)  =  A(n\  Defining  the  error  generator 
matrix,  D(Ar',  such  as 

DW=I-H(1V>  (10) 

EW  is  directly  obtained  from  A:  AD(Ab  =  EW. 

(2)  is  presented  as  the  QR  factorization,  e.g.,  [12],  in 
which  (5)  does  not  generally  take  place,  either  in  the 
GSO,  but  5=  (1,2,...  ,7). 

As  the  evaluation  function,  a  normalized  version  of 
the  squared  Euclidean  norm  of  is  employed  and 

denoted  by 


A  is  reconstructed  (or  approximated)  from  N 
columns  selected  from  A.  The  reconstructed  matrix  is 
expressed  as  A{N )  where  the  superscript  ( N )  indicates 
the  number  of  selected  columns  in  A.  The  reconstruc¬ 
tion  error  matrix  E  is  defined  as  E*^  =  A  —  A^ 
and  it  consists  of  error  vectors  e*= \j: 


Ew 


(4) 


N  also  denotes  the  iteration  number.  The  iterative 
process  selects  one  column  in  A  in  each  iteration.  As 
the  initialization  (when  N  =  0),  E^0)  is  set  equal  to  A. 

In  each  iteration,  N  =  1  :  I,  the  index  of  error  vec¬ 
tor  whose  2-norm  is  the  largest  of  is  detected 

and  stored  as  Nth  element,  sn  ,  in  the  lead  selection 
sequence  list  S: 

sjv  =  argmaxHe^-1^^  ,  N  =  1:  I  (5) 

i—1 

S  =  (si,  s2,  ...,  si)  (6) 

qN  is  set  equal  to  detected  by  (5): 

qN  =  ,  N  =1:1  (7) 


and  row  N  in  R  are  set  by 


VN,i  = 


t  (JV-l) 

<In te) 

e(jv-i)te(jv-i)  ’ 


i  =  1  :  I 


(8) 


e(N)  = 


l|E(JV)||g 

l!AIU 


(ii) 


e^N'  is  the  “overall”  reconstruction  error  ratio;  i.e.,  it 
does  not  indicate  an  error  level  of  particular  lead. 

HW  and  DW  can  be  applied  to  any  matrix  if  its 
matrix  size  is  the  same  as  A  from  which  the 
and  the  D(Ar)  are  computed.  The  process  to  compute 
them  is  referred  to  as  training ,  and  evaluation  by 
applying  to  a  matrix,  which  may  or  may  not  be 

the  A,  is  referred  to  as  test  in  the  following. 

Four  tests  were  executed  in  each  model  such  that 
e(N) )  tv  =  1  :  J,  were  computed  applying  the 

D i'n\N  =  1:1,  obtained  from  all  the  four  models. 


IV.  Result 

el'N')  in  the  tests  with  different  models  are  compared 
in  Figure  1.  The  number  of  leads,  N,  was  incre¬ 
mented  from  1  to  120  where  120,  which  is  the  num¬ 
ber  of  full-set  BSM  leads.  The  minimum  N  to  obtain 
e(N'  less  than  or  equal  to  certain  values  are  briefly 
summarized  in  Table  1.  They  are  not  significantly  dif¬ 
ferent  from  the  numbers  previously  reported  by  other 
researchers  [6,  7,  9].  The  cardiac  phase  of  the  heart 
did  not  strongly  affect  the  test  performance;  i.e.,  the 
same  electrode  set  and  reconstruction  coefficients  may 
be  used  in  models  of  different  cardiac  phases  if  all  the 
other  factors  in  the  models  are  same.  The  inhomo¬ 
geneities  made  significant  differences,  especially  when 


min  N  for  efN')  <  0.1 

Test  Model 

Training  Model 

i-sys 

i-dia 

h-sys 

h-dia 

i-sys 

16 

18 

15 

18 

i-dia 

17 

17 

15 

16 

h-sys 

18 

20 

12 

14 

h-dia 

20 

21 

13 

13 

min  N  for  <  0.01 

Test  Model 

Training  Model 

i-sys 

i-dia 

h-sys 

h-dia 

i-sys 

30 

35 

37 

37 

i-dia 

31 

32 

35 

38 

h-sys 

53 

48 

29 

31 

h-dia 

45 

45 

29 

30 

Table  1:  The  minimum  numbers  of  leads,  TV,  to  obtain  recon¬ 
struction  error  ratio  >  less  than  or  equal  to  0.1  and  0.01  are 
tabulated.  Lead  selection  sequences  5  and  lead  reconstruction 
error  generator  matrices  D '  -''v *  were  obtianed  in  traninigs.  In 
the  tests  were  computed  using  the  D 1 N 1  for  TV  =  1.  3, . . . 
as  simulations  of  incrementing  the  leads  according  to  S. 

e1' ,v-)  became  low.  The  difference  of  conductivities  ex¬ 
amined  here  is  an  extreme  case  —  inhomogeneous  or 
infinite  extent  of  homogeneous  conductivity.  The  re¬ 
sults  suggest  that  the  inhomogeneity/homogeneity  set¬ 
ting  should  not  be  changed  for  accurate  lead  recon¬ 
structions  in  model-based  studies. 

The  reconstruction  matrix,  H,  does  not  have  an  ef¬ 
fect  of  amplification  or  attenuation.  Its  effect  is  a  rota¬ 
tion  in  a  vector  space.  Thus,  only  relative  conductivi¬ 
ties  in  each  model  concerns  when  models  are  different. 
The  effects  by  changing  models  are  most  likely  due  to 
distribution  patterns  of  conductivities. 

V.  Conclusion 

The  numbers  of  leads  necessary  for  certain  levels  of 
full-BSM  reconstruction  were  not  significantly  differ¬ 
ent  from  those  reported  by  other  researchers.  Differ¬ 
ence  of  cardiac  phases  of  the  models  did  not  affect  the 
reconstruction  performance  while  difference  of  inhomo¬ 
geneous/homogeneous  conductivities  clearly  affected. 
The  lead  selection  method  in  the  preset  study  is  math¬ 
ematically  the  same  as  one  applied  to  real  ECG  dataset 

[7].  The  next  interest  will  be  exchangeability  of  lead 
selection  and  reconstruction  coefficients  obtained  by 
model  approach  and  by  dataset  approach. 
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Figure  1:  The  error  generator  matrices,  were  obtained  from  four  thorax  models  (as  training  models)  and  applied  to  the  four 

models  (as  test  models).  The  overall  error  ratio  is  an  evaluation  function  of  lead  reconstruction  performance.  N  is  the  number  of 
leads  to  be  selected,  i.e.,  it  is  the  number  of  vectors  selected  from  the  transfer  matrix  of  test  model  from  which  the  entire  transfer 
matrix  itself  is  reconstructed  using  the  reconstruction  matrix.  The  lead  selections  and  incorporating  reconstruction  matrices  were 
obtained  from  training  models.  The  error  generator  matrix  generates  the  error  matrix  in  reconstruction  tests.  The  plots  above  show 
changes  in  the  reconstruction  error  by  incrementing  the  number  of  selected  leads.  The  abbreviations  of  the  models  are  as  follows: 
i-sys:  inhomogeneous  end-systolic  model,  i-dia:  inhomogeneous  end-diastolic  model,  h-sys:  homogeneous  end-systolic  model,  h-dia: 
homogeneous  end-diastolic  model. 


